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Summary 


Bat species around the world have recently been recognized as major reservoirs of 
several zoonotic viruses, such as severe acute respiratory syndrome coronavirus 
(SARS-CoV), Middle East respiratory syndrome coronavirus (MERS-CoV), 
Nipah virus and Hendra virus. In this study, consensus primer-based reverse tran- 
scriptase polymerase chain reactions (RT-PCRs) and high-throughput sequencing 
were performed to investigate viruses in bat faecal samples collected at 11 natural 
bat habitat sites from July to December 2015 in Korea. Diverse coronaviruses were 
first detected in Korean bat faeces, including alphacoronaviruses, SARS-CoV-like 
and MERS-CoV-like betacoronaviruses. In addition, we identified a novel bat 
rotavirus belonging to group H rotavirus which has only been described in 
human and pigs until now. Therefore, our results suggest the need for continuing 
surveillance and additional virological studies in domestic bat. 


et al., 2014; Yang et al., 2014; Mortlock et al., 2015). Reas- 
sortant group A rotavirus was also detected in the faeces of 


Bats are considered reservoirs of several emerging infec- 
tious disease. Recently emerging human viruses, such as 
severe acute respiratory syndrome coronavirus (SARS- 
CoV), Middle East respiratory syndrome coronavirus 
(MERS-CoV), Nipah virus, Hendra virus and Ebola virus, 
are thought be bat-borne viruses (Han et al., 2015). In 
addition, new species of influenza A viruses, lyssaviruses, 
paramyxoviruses and coronaviruses continue to be discov- 
ered in bats around the world (Tong et al., 2013; Banyard 
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a straw-coloured fruit bat (Esona et al., 2010). Therefore, 
viruses that originate in bats may be a source of additional 
spill over from wildlife into domestic animals and humans 
(Plowright et al., 2015). 

Since the outbreak of MERS-CoV in Korea in 2015, the 
prediction and control of newly emerging viruses has 
gained urgency. However, in Korea, there have been a few 
studies on viruses in domestic bat species. Serological evi- 
dence of Japanese encephalitis virus (JEV) was reported in 
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the early 1990s (Lee and Lee, 1992). In 1995, there was a 
report on the isolation and genetic characterization on a 
hantavirus from bats (Jung and Kim, 1995). More informa- 
tion of viruses in bats is needed to clarify the role of these 
animals in infectious diseases epidemiology. 

Therefore, in this study, we investigated viruses in bat 
species in Korea, using 49 faecal samples collected from July 
to December 2015 in 11 sites in natural bat habitats. 


Materials and Methods 


Samples 

From July to December in 2015, 49 faecal samples were col- 
lected at 11 sites in natural bat habitats, such as caves, an 
abandoned mine and under a bridge (Table 1). Fresh bat 
faeces were placed into transport medium in a 10% suspen- 
sion and were transported to the laboratory for further lab- 
oratory analysis. 

The major bat species at the collection sites were deter- 
mined based on morphology and on previous data from 
bats’ roosting sites (Han et al., 2012). In this study, the 
major bat species in the collecting sites were Rhinolophus 
ferrumequinum, Myotis macrodactylus, Myotis aurascens 
Kuzyakin, Myotis petax and Miniopterus schreibersii. 
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Viral detection 


Consensus primer-based reverse transcriptase polymerase 
chain reactions (RT-PCRs) were performed to detect 
influenza A viruses, coronaviruses, lyssaviruses and param- 
yxoviruses (Fouchier et al., 2000; Poon et al., 2005; 
Vazquez-Moron et al., 2006; Tong et al., 2008). The 
detailed information on the RT-PCRs is in Table S1. A total 
of 49 bat faecal samples were tested. The amplicons from 
positive RT-PCRs were sequenced using target-specific for- 
ward and reverse primers synthesized by Cosmogenetech 
Co. Ltd (Daejeon, Korea). The nucleotide data obtained 
from Sanger sequencing of the PCR fragments were further 
analysed with related sequences in GenBank using BioEdit 
(Hall, 1999) and MEGA version 6 (Tamura et al., 2013). 


High-throughput sequencing 


In total, 16 bat faecal samples that were collected in four 
sites were used for high-throughput sequencing (Table 1). 
The supernatants from suspensions of one, four, five and 
six faecal samples collected in Danyang, Wonju, Taebaek 
and Uljin, respectively, were pooled by sampling site. Each 
of the four pooled samples was filtered through a 0.2-um 


Table 1. Information on bat faecal samples 
Cases Collected date Region Site name Major bat species No. of samples 
1 24 July 2015 Danyang Under HSA Bridge Unknown 6° 
2 4 August 2015 Wonju HO cave Rhinolophus ferrumequinum, 4c 
Myotis macrodactylus 
3 4 August 2015 Taebaek YY cave* Myotis aurascens Kuzyakin 5° 
Myotis petax 
Rhinolophus ferrumequinum 
4 6 August 2015 Uljin SL cave? Unknown 1° 
5 9 October 2015 Jeongsun PG cave Rhinolophus ferrumequinum 4 
6 4 October 2015 Inje BM abandoned Rhinolophus ferrumequinum 6 
mine 
7 18 October 2015 Taebaek YY cave Myotis aurascens Kuzyakin 3 
Myotis petax 
Rhinolophus ferrumequinum 
8 25 October 2015 Bonghwa DRN cave Rhinolophus ferrumequinum 3 
Myotis macrodactylus 
9 18 November 2015 Taebaek YY cave Myotis aurascens Kuzyakin il 
Myotis petax 
Rhinolophus ferrumequinum 
10 19 November 2015 Uljin SL cave? Unknown 1 
11 1 December 2015 Hapcheon BT cave Unknown 3 
12 1 December 2015 Munkyung BGS cave Miniopterus schreibersii 4 
13 December 2015 Yeongwol YJ cave Rhinolophus ferrumequinum 4 
14 December 2015 Miwon CS cave Rhinolophus ferrumequinum 4 


%collected at the same site (YY cave) on a different date. 
®collected at the same site (SL cave) on a different date. 
“pooled and sequenced by high-throughput sequencing. 
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filter and ultracentrifuged at 252 000 g for 1 h. Each pellet 
was resuspended in 500 pl of 1x digestion buffer (Turbo 
DNA Free Kit; Ambion, Darmstadt, Germany), and a 
total of 10 yl of Turbo DNase was added. After incubat- 
ing for 30 min at 37°C, suspensions were transferred to 
a 1.5-ml reaction tube and centrifuged at 1500 g for 
3 min. The supernatants were used for RNA extraction. 
RNA was extracted using Trizol LS (Invitrogen Corp., 
Carlsbad, CA, USA), following the manufacturer’s man- 
ual. The extracted RNA was submitted to Macrogen 
(Seoul, Korea) for high-throughput sequencing in a 
HiSeq 2000 sequencing system based on the transcrip- 
tome de novo sequencing platform. 

The obtained viral contigs were analysed and annotated 
by the MG-RAST server (Meyer et al., 2008). The cut-off 
for the annotation was 10~°, 60% and 15 for maximum 
e-value, minimum percentage identity and minimum align- 
ment length, respectively. The contigs annotated as mam- 
malian viruses were validated through BLASTn (http:// 
www.ncbi.nlm.nih.gov). The GenBank accession numbers 
for the nucleotide sequences in this study are KU528584— 
KU528594. 


Results 


PCR-based detection of influenza A viruses, 
coronaviruses, lyssaviruses and paramyxoviruses in bat 
faeces 

There were no samples positive for influenza A viruses, lys- 
saviruses or paramyxoviruses. 

Among 14 cases of the submitted samples, three cases 
were positive for coronavirus by RT-PCR, showing the 
expected 440 nt target. These included the following: (Case 
no. 1) one of four samples from HO cave in Wonju where 
Rhinolophus ferrumequinum and Myotis macrodactylus were 
the major species present, (Case no. 6) one of six samples 
BM abandoned mine in Inje where Rhinolophus ferrume- 
quinum was the major species found and (Case no. 12) two 
of four samples BGS (Case no. 12) cave in Munkyung 
where Miniopterus schreibersii was found to be the major 
species (Table 1). 

In the HO cave, one of four faecal samples (B15-8) was 
positive by coronavirus-specific RT-PCR. When the PCR- 
positive fragment was sequenced, it was found to be a bat 
coronavirus (Bat CoV) belonging to the genus alphacoron- 
avirus and showing 89% nucleotide sequence identity with 
Bat CoV A620 (GenBank accession no. DQ648828) 
detected in China in 2005. 

In the BM abandoned mine, one of four faecal samples 
(B15-21) was found to be a Bat CoV belonging to the genus 
betacoronavirus and showing 98% nucleotide identity with 
SARS-like Bat CoV Lingbao-331 (GenBank accession no. 
KF294456) detected in China in 2012. 
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In the BGS cave, two of four faecal samples (B15-40 and 
B15-41) were found to be Bat CoVs belonging to the genus 
alphacoronavirus and showing 99% nucleotide identity 
with Bat CoV A895 (GenBank accession no. DQ648841) 
and 98% nucleotide identity with Bat CoV A910 (GenBank 
accession no. DQ648849), both of which were detected in 
China in 2005. 


High-throughput sequencing from the pooled bat faecal 
samples 


As shown in Table 1, cases 1, 2, 3 and 4 were further anal- 
ysed by high-throughput sequencing. By Hiseq 2000 
sequencing, 4341, 162 559, 7690 and 16 657 contigs were 
obtained from cases 1, 2, 3 and 4, respectively. Based on the 
annotation by the MG-RAST system, 
sequences that are related to mammalian viruses, insect 
viruses, plant viruses, fungal viruses and bacteriophages 
were identified. 

The mammalian viral sequences annotated by the MG- 
RAST system were further validated by BLASTn (http:// 
www.ncbi.nlm.nih.gov). The BLASTn result is presented 
in Table 2. Fifteen viral sequences from case no. 1 were 
similar to porcine rotavirus SKA-1 (69-73% nucleotide 
identity), Bat CoV SC2013 (70-91% nucleotide identity), 
Bat CoV HKU5-5 (82% nucleotide identity), Alphacoron- 
avirus Eptesicus fuscus-related strains (95-96% nucleotide 
identity) and Rhinolophus pearsoni bunyavirus Shaanxi 
2011 (72% nucleotide identity). Three viral sequences 
from case no. 2 were similar to Banna virus 02VN078b 
(90% nucleotide identity), Bat CoV CDPHEI15 (80% 
nucleotide identity) and Bat CoV Neixiang-31 (82% 
nucleotide identity). There was no mammalian viral 
sequences in cases 3 and 4. 


several viral 


Genetic characteristics of bat coronaviruses in Korea 


Partial RNA-dependent RNA _ polymerase (RDRP) 
sequences obtained from B15-8 (sample in case no. 2), 
B15-21 (sample in case no. 6), B15-40 and B15-41 (sample 
in case no. 12) were phylogenetically compared with other 
reference sequences (Fig. la). Three of these sequences 
clustered with the alphacoronaviruses (B15-8, 40 and 41), 
and the other one clustered with the betacoronaviruses 
(B15-21). Bat CoV B15-8 grouped with bat CoVs found in 
China, showing 85.7-88.8% nucleotide identities, and was 
found to be closely related with Bat CoV-512 and porcine 
epidemic diarrhoea virus (PEDV)-CV777, showing 78.6% 
and 74.6% nucleotide identities, respectively. Bat CoV B15- 
21 grouped with SARS-CoV-like Bat CoVs, such as Bat 
CoV Rfl and Bat CoV HKU3 with 90.3-97.7% nucleotide 
identities, and was found to be closely related to human 
SARS-CoVs (SARS-CoV Tor2, MA15, ExoN1 strains), 
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BLASTn search parameters 


Matched strain 
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Table 2. Mammalian viral sequences matched 
with contigs from high-throughput sequencing 


Porcine rotavirus SKA-1 
Porcine rotavirus SKA-1 
Porcine rotavirus SKA-1 
Porcine rotavirus SKA-1 
Porcine rotavirus SKA-1 
Porcine rotavirus SKA-1 
BtVs-BetaCoV/SC2013 
BtVs-BetaCoV/SC2013 
Bat coronavirus HKU5-5 
BtVs-BetaCoV/SC2013 
Alphacoronavirus Eptesicus 


fuscus/Appalachian Ridge/ 
P1-C1148/IT/USA/2009 


Alphacoronavirus Eptesicus fuscus 


strain NeCoV isolate Ffrag_ 
SAGA_bigbrown_2010 


Alphacoronavirus Eptesicus 


fuscus/Appalachian Ridge/ 
P1-C837/IT/USA/2009 


Anlong Mss bat coronavirus 


isolate Neixiang-32 


Rhinolophus pearsoni bunyavirus 


isolate Bat Ro-BunyaV/ 


Banna virus isolate O2VNO78b 
Bat coronavirus CDPHE15/USA/ 


Size % Nucleotide 

Contigs* (bases) F£-value identity Accession No. Strain name 
C6091 404 2e-37 73 AB576634 

C6719 438 Je-21 69 AB576631 

C7554 211 9e-19 73 AB576625 

C9980 215 5e-22 73 AB576629 

c10008 = 285 2e-35 73 AB576629 

C10298 = 355 2e-48 73 AB576631 

C3339 230 Je-55 83 KJ473821 

C6171 203 1e-16 70 KJ473821 

c10151 219 1e-36 82 EFO65512 

C10579 =. 207 3Be-67 991 KJ473821 

C6073 243 6e-98 95 HQ585081 

C8820 417 0.0 96 JX537911 

C8271 270 2e-61 95 HQ585082 

C7729 209 2e-21 76 KF294382 

C5770 313 2e-34 72 KC154063 

Shaanxi2011 

C700585 207 8e-70 90 EU265697 

C817114 201 6e-40 80 KF430219 

2006 

C695165 226 6e-53 82 KF294381 


Neixiang Md bat coronavirus 


isolate Neixiang-31 


*Contigs obtained from case no. 1 and case no 2. No mammalian viral sequences were detected in 


case no. 3 and case no 4. 


showing 89.1% nucleotide identity. Bat CoV B15-40 and 
B15-41 grouped with Bat CoV 1A and 1B, and nucleotide 
identities within the group were 92.3-97.1%. 

Because an RDRP sequence from _betacoronavirus- 
related Bat CoV B15-1-6 (detected from the pooled sample 
in case no. 1) was not obtained, a partial spike gene 
sequence (c10151) was compared with the spike genes of 
reference betacoronaviruses (Fig. 1b). In addition, the full 
spike gene sequence obtained from Bat CoV B15-21 was 
also analysed. Bat CoV B15-1-6 grouped with MERS-CoV- 
like Bat CoVs (Bat CoV strains HKU5, A434, SC2013 and 
HKU4), showing 80.9-82.2% nucleotide identities, and was 
found to be closely related to MERS-CoV EMC, showing 
76.6% of nucleotide identity. 

The amino acids sequence encoded by the full spike gene 
of B15-21 was aligned and compared with human SARS- 
CoV and SARS-CoV-like Bat CoV sequences (Figure S1). 
The S2 domain is more conserved than the $1 domain in 


the spike protein. The receptor binding domain (RBD) of 
the SARS-CoV-like Bat CoVs, except for RaSHC014, had 
two major deletion sites, TGNYN (446-450) and 
PFSPDGKPCTPPA (472-484), compared to human SARS 
coronaviruses sequences. The receptor binding motif 
(RBM) and human ACE2 (hACE2)-interacting residue 
were highly variable between human SARS-CoVs and 
SARS-CoV-like Bat CoVs but RaSHC014 strain. 


Phylogenetic characterization of Group H bat rotavirus in 
Korea 


Porcine rotavirus-related sequences from pooled bat faeces 
in case no. 1 were analysed with reference sequences of 
rotaviruses in human and animals (Fig. 2). Based on the 
partial sequences of the VP1, VP3 and VP4 regions, the 
rotavirus-related sequences were found to belong to rota- 
virus group H. The group H bat rotavirus in this study 
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40 -— BatCoVv-A620(00648828) 
(a) 45[- | Batcov-a613(00648826) (b) 
53 BatCoV-HKU6(D0249224) 

fico BatCoV-A6897(DQ648842) 

BatCoV-A701(D0648833) 

[— BatCoV-A632(00648830) 

‘99 |_— BatCov-A633(D0648831) 

158 

BatCoV-512(NC_009657) 

PEDV(KC189944) 

14) (_________gatcov-aju21(kc779226) 

[— HCoV-NL63(D0445911) 

1 L_— HCoV-229E(JX503061) 

[— BalCoV-RoHKU10(NC_018871) 

7100 BatCov-HiHKU10(J0989273) 

BatCoV-1B(NC_010436) 

77 00 B15-41 
B15-40 

51 L— BatCoV-1A(NC_010437) 

BatCoV-HKU2(NC_009988) 

FCoV(EU186072) 

(— TGEV(AJ271965) 

1001 ccow(kP9si644) 

400 -— !BVINC_001451) 

L_ Tev(nc_010800) 

70 -— BCV(U00735) 

100 L__PHEWINC_007732) 

100 HCoV-0C43(NC_005147) 

09] [— HCoV-HKU1(NC_008577) 

71 MHV(AB551247) 

BatCoV-HKUS(NC_009021) 

100 -— MERS-CoV(NC_019843) 


100 


Alpha CoVs 


Gamma CoVs 
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Bat CoV-HKU5(KC522101) 
Bat CoV-A434(00648790) 
Bat CoV-SC2013(KJ473821) 


MERS CoV-like 
Clade 
MERS-CoV(NC_019843) 

Bat CoV-HKUS-1(NC009019) 
Bat CoV-RsSHC014 (KC88 1005) 
SARS CoV-Sino1-11(AY485277) 
SARS CoV-Urbani (AY278741) SARS CoV-like 
————— Bat CoV-Rp3 (00071615) Clade 

ro) 615-21 

7 Bat CoV-RI(004 12042) 
Bat CoV-A701(00648793) 
HCoV-NL63(AY567487) 


Beta CoVs 


| Batcov-HkU4(NC_009019) 
35 -— SARS CoV-Tor2 (JX163928) 

100[  |__ SARS Cov-Ma15 (JF292915) 

SARS CoV-ExoN1 (FJ882955) 

4100 -— BatCov-HkU3(00022305) 

|__ SARS-like Bat CoV-HKU3-1 (00022305) 

a Bi5-21 

Bat CoV-Rp (JX993987) 

Bat CoV-273 (00648856) 

[— SARS-like Bat CoV-Lingbao-331 (KF294456) 

‘92 L__ SARS-lke Bat CoV-RA1 (00412042) 


BeTA CoVs 


700 


100 


HCoV-229E (AF 304460) 
PEDV-CV777(AF 353511) 
BatCoV-A5 15(0Q648791) 
BatCoV-HKU10(JQ989270) 


Alpha CoVs 


Fig. 1. Phylogenetic trees constructed from bat coronavirus sequences detected in this study and from other coronaviruses. (a) bootstrap consensus 
tree inferred from 1000 replicates with the UPGMA method (evolutionary distance parameter: maximum composite likelihood method) based on 351 
bases of the RNA-dependent RNA polymerase gene, (b) bootstrap consensus tree inferred from 1000 replicates using the neighbour-joining method 
(evolutionary distance parameter: maximum composite likelihood method) based on 163 bases of the spike gene. The box indicates the bat CoV 


sequences obtained in this study. 


showed 69.4—72.7%, 71.1-72.9% and 67.3—72.2% nucleo- 
tide identities with the other group H rotaviruses, based on 
VP1, VP3 and VP4 partial sequences, respectively. 


Discussion 


According to a recent report, there are 23 species of bats in 
Korea (Han et al., 2012). Most Korean bats are thought to 
be insectivores, as no fruit bats have been found in Korea 
to date, except for some fruit bats as pets or indoor exhibi- 
tions (Han et al., 2010). Through the target-specific RT- 
PCR, partial RDRP sequences of three alphacoronaviruses 
(B15-8, 40 and 41) and one betacoronavirus (B15-21) were 
detected in the faeces of Korean bats. Bat CoV B15-8 was 
detected in faeces collected at the HO cave where Rhinolo- 
phus ferrumequinum and Myotis macrodactylus were the 
major species. This virus was found to be closely related 
with Bat CoVs recently found in China and PEDV, a por- 
cine diarrhoeic virus around world. On the other hand, Bat 
CoVs B15-40 and 41 were detected at the BGS cave, where 
Miniopterus schreibersii was the major species. These 
sequences were closely related with Bat CoV 1A and 1B, 
which were reported to have a clear host species restriction 
to Miniopterus magnate and Miniopterus pusillus in Hong 
Kong (Chu et al., 2008). As Bat CoVs B15-40 and 41 were 
also found in faeces in a cave where Miniopterus schreibersii 
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was dominant, Bat CoV strains in the same lineage with 
Bat CoV 1A and 1B may be adapted to the Miniopterus spp. 

It was interesting that Bat CoV B15-21 found in this 
study was closely related to SARS-CoV-like Bat CoVs, such 
as Bat CoV Rfl, based on the partial RDRP and S genes, 
respectively. The newly identified Bat CoV B15-21 clustered 
with the SARS-CoV-like Bat CoV group in the phylogenetic 
analysis. The major bat species in the site where Bat CoV 
B15-21 was collected was Rhinolophus ferrumequinum, 
which is consistent with a previous report that SARS-CoV- 
like bat CoVs were mostly found in Rhinolophus spp. (Yuan 
et al., 2010). The amino acids of its full spike gene showed 
that the B15-21 strain is closely related to SARS-CoV-like 
Bat CoVs Rp3, Rfl, HKU3-2 and 273 rather than human 
SARS-CoVs and Bat CoV RaSHC014, which was reported 
to have high potential for human emergence (Menachery 
et al., 2015) and showed two specific deletion sites in the 
RBD of the spike protein. One additional betacoronavirus 
(B15-1-6), which was not detected by the consensus pri- 
mer-based RDRP-specific RT-PCR (Poon et al., 2005), was 
detected by high-throughput sequencing as shown in 
Table 2. Three contigs were similar to Bat CoV SC2013, 
which is a MERS-related betacoronavirus that was found in 
bats in China (Yang et al., 2014). One contig was also simi- 
lar to Bat CoV HKU5-5, which clustered with MERS-CoV 
(Woo et al., 2007; Annan et al., 2013). Phylogenetic 
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(a) 42, Bat RVH (c10008} 
51 ‘SKA-1(Porcine)-AB576629 
Pl J19(Human)-DQ113897 RVH 
90 B219(Human)-EF453355 
03V0567(Chicken)-NC_021590 RVG 
CAL-1(Human)-EU490414 
700 Bang373(Human)-NC_021541 RYE 
43 03V0568(Chicken)-NC_021625 
05V0049(Chicken)-NC_014511 RVD/RVE 
00 400 KU(Human)-AB022765 RVA 
OSU(Porcine)-GU199514 
36| Bristol(Human)-NC_007547 RVC 
Er) Cowden(Porcine)-M74216 
(b) 100; J19(Human)-DQ113900 
57 B219(Human)-EF453357 
95 |__ska.-1(Porcine)-AB576631 RVH 
Bat RVH (C10298)| 
——___080567(Chicken)-NC_021589 RVG 
at Cal-1(Human)-EU490417 RVB 
Tod Bang-373(Human)-NC_021551 
400 KU(Human)-AB022767 
OSU(Porcine)-AY277921 RVA 
700 99 Bristol(Human)-NC_007574 RVC 
Cowden(Porcine)-M74219 
32 05V0049(Chicken)-NC_014514 RVD/RVE 
58 (03V0568(Chicken)-NC_021630 
100 J19(Human) 
(c) 84 B219(Human) 
BR59(Porcine) 
99 SKA-1(Porcine) 
03V0567(Chicken) RVG 
5 IDIR(Murine) 
aa -———— ADRV(Human) RVB 
700 Cal-1(Human) 
99 Bang-373(Human) 
4 KU(Human) 
OSU(Porcine) ne 
wad 84 Bristol(Human) RVC 
Cowden(Porcine) 
70 05V0049(Chicken) RVD/RVF 
58 


03V0568(Chicken) 


Fig. 2. Phylogenetic trees constructed from Group H rotavirus-related 
sequences detected in bat faeces in this study and sequences of other, 
representative rotavirus strains. The bootstrap consensus tree inferred 
from 1000 replicates using the neighbour-joining method (evolutionary 
distances were computed using the Kimura 2-parameter method). (a) 
275 bases of the VP1 region, (b) 347 bases of the VP3 region, (c) 202 
bases of the VP4 region. The box indicates the group H bat rotavirus 
sequences obtained in this study. 


analysis based on partial spike gene sequences of Bat CoV 
B15-1-6 showed that it belonged to a MERS-CoV-like clade 
with previously known MERS-CoV and MERS-CoV-like 
Bat CoVs. 

Although information of the associated bat species was 
not clearly available, this study first report that SARS-CoV 
and MERS-CoV-like bat CoVs exist in Korea. This may 
imply that other strains or types of those coronaviruses 
may circulate in bat species in Korea. Therefore, we cannot 
ignore appearance of novel or variant coronavirus through 
local emergence in Korea, like SARS and MERS in China 
and Middle East, respectively. 

In addition to coronaviruses, this study is the first to 
determine the existence of group H rotaviruses in bat faeces 
collected in Korea. Group H rotavirus is a recently pro- 
posed group of rotaviruses that include strains ADRV-N 
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and B219, which infect human adults (Alam et al., 2007; 
Matthijnssens et al., 2012). So far, group H rotaviruses 
have only been reported in human and pigs (Molinari 
et al., 2015), but this study provides evidence that bat spe- 
cies may be a host of group H RVs. To confirm that, there 
should be follow-up studies including virus isolation and 
characterization, genomic analysis, continuous surveillance 
and VP6-based classification (Matthijnssens et al., 2012) to 
find its prevalence, epidemiology and zoonotic potential. 

Localized emergences of zoonotic diseases are usually 
due to expansion of the wildlife-human interface (Morse 
et al., 2012). The human—bat interface in Korea has not 
been fully described. People usually recognize that bats are 
flying and roosting in the deep forest or caves. However, 
the sites where the bat faecal samples were collected showed 
some evidence of human-—bat contact (Figure S2). Some 
caves or abandoned mine where bats are roosting are 
located near human habitats. People living there frequently 
visit into the caves or the abandoned mine to take a rest or 
to ferment foods. Some caves were also used for some reli- 
gious ceremony by Korean traditional shamans. 

In this study, SARS-CoV-like and MERS-CoV-like bat 
CoVs and group H rotavirus were detected for this first 
time in Korea, which may be of interest because of their 
zoonosis potential. Several kinds of evidence for human— 
bat contact were also found. Therefore, the potential for 
novel or variant human viruses from bats in Korea should 
not be ignored. Continuous surveillance and additional 
virological research should be a priority to prevent future 
emergences of zoonotic viral diseases. 
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